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Abstract 

Angular correlation of jets produced in association with a massive scalar, vector or 
tensor boson is crucial in the determination of their spin and CP properties. We study 
jet angular correlations in events with a high mass bottom quark pair or a top quark pair 
and two jets at the LHC, whose cross-section is dominated by the virtual gluon fusion 
sub-processes when appropriate kinematic selection cuts (vector-boson fusion cuts) are 
applied. We evaluate helicity amplitudes for sub-processes initiated by qq, qg and gg 
collisions in the limit where the intermediate gluons are collinear to the initial partons. 
We first obtain a general expression for the azimuthal angle correlations among the 
dijets and ti or 66, in terms of the gg — )• ti or 66 helicity amplitudes in the real gluon 
approximation of the full matrix elements, and find simple analytic expressions in the 
two kinematic limits, the production of a heavy quark pair near the threshold, and in 
the relativistic limit where the invariant mass of the heavy quark pair is much larger 
than the quark mass. For 66 + 2 jets we find strong azimuthal angle correlations which 
are distinct from those expected for events with a CP-even or odd scalar boson which 
may decay into a 66 pair. For ti -\- 2 jets we find that the angular correlations are 
similar to that of a CP-odd scalar+2 jets near the threshold Mjf ~ 2?TT,t, while in the 
relativistic limit they resemble the distribution for 66 + 2 jets. These correlations in 
the standard QCD processes will help establish the experimental technique to measure 
the spin and CP properties of new particles produced via gluon fusion at the LHC. 
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1 Introduction 



Angular correlations among jets produced in association with scalar bosons and massive 
gravitons have been the subject of several studies [HE]. It has been found that experimental 
determination of such correlations can give us important information about the spin and 
CP-properties of these particles. For example, by studying the distribution of azimuthal 
angle difference between two tagging jets (A0jj) in the gluon fusion (GF) production of a 
scalar (H) or pseudo-scalar (A) Higgs with two jets, it has been shown that the CP-odd 
and CP-even Higgs can clearly be discriminated at the LHC. It has also been shown that 
for jets produced with massive gravitons (G), the distribution is flat while the sum 

of the azimuthal angles show a characteristic distribution |2]. It was further demonstrated 
analytically in Ref. |2] that these angular correlations arise from the quantum interference 
of different helicity states of the intermediate vector bosons, which can either be a weak 
boson or a gluon. Another critical observation made in Ref. [2] is that the amplitudes for 
H/A/G + 2 jets production can be approximated by those of the vector boson fusion (VBF) 
sub-processes when a specific set of kinematic cuts are applied, where the two tagging jets 
are demanded to be in the opposite hemispheres of the detector and also to be well separated 
in rapidity (the so called VBF cuts), and in addition their transverse momenta are subjected 
to a slicing cut, restricting them to be sufficiently low as compared to the H/A/G masses. 

In this paper, we study heavy quark pair QQ {Q = b, t) production at the LHC in 
association with two jets, and study the azimuthal angle correlations among these jets. 
Our goal is to predict the jet angular correlations for the SM processes by using the same 
approximation, so that the experimental technique to measure the new particle properties 
via initial state radiation patterns can be established using the ample SM processes. We 
first calculate the helicity amplitudes for gluon fusion processes initiated by gg, qg, gg initial 
states at the LHC leading to QQ + 2 jets production. Here, we work in the limit in which the 
intermediate gluons are collinear to the initial partons. We work out the angular correlations 
using the formalism of Ref. [2] , and then verify our approximate analytical results in several 
kinematic regions by comparing them with the exact parton level matrix elements in the 
tree level. 

This paper is organized as follows. In section 2 we briefly review the helicity formalism 
used to obtain our analytical results. In section 3 we give the spin summed matrix-element 
squared for QQ + 2 jets production via gluon fusion in the collinear limit for the exchanged 
gluons, and then analyze the results in two different kinematic limits, namely that of the ti 
production threshold [Q = t) in which the invariant mass of the ti pair, Mjj, is very close to 
2mt, and that of the relativistic limit where Mqq is much higher than 2mQ, for Q = b and 
t. In section 4 we show the exact numerical results for these angular correlations at parton 
level, where the expected azimuthal angle correlations are reproduced once appropriate final 
state selection cuts are applied. Finally, we summarize our findings in section 5. 

2 Helicity formalism : brief review 

The helicity amplitude formalism for vector boson fusion processes has been discussed in 
detail in Ref. [2] in the context of production and decay of a heavy particle in association 
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Figure 1: Schematic Feynman diagram for the parton level QQjj production process in 
the on-shell gluon approximation, where denote quarks, anti-quarks or gluons. The four- 
momenta of each particle are also shown along the particle lines. 



with two jets. We follow the same notations and conventions in this paper. The inclusive 
process. 



can proceed, via gluon fusion, through the sub-processes aia2 — ?■ a3a4g*g* — )■ 0304(5(5, where 
aia2 is one of the three possible initial states {qq, qg, gg}, g* is a t-channel intermediate gluon 
and q stands for a quark or anti-quark of any flavour: 



In addition to gluon fusion, the above sub-processes receive contributions from all other 
diagrams at the same order in perturbation theory to make a gauge-invariant physical ampli- 
tude. However, in this and the next section, we shall consider only the gluon fusion diagrams. 
As mentioned earlier, it has been demonstrated in Ref. [2j that after applying the VBF se- 
lection cuts, the GF contribution can be made to dominate the total cross-section. We shall 
also show the validity of this approach by first determining analytically the distributions 
predicted by the GF diagrams in various kinematic limits, and then comparing them with 
the exact matrix elements. 

To begin with, we define a common set of kinematic variables for the QQjj sub-processes 
in Eqn. |2] as 



The helicities of the quarks or anti-quarks have the value cTj/2 while for on-shell gluons they 



PP ^ QQjj + anything , 



(1) 



qq qqQQ 

Qg QgQQ 

99 99QQ 



(2a) 
(2b) 
(2c) 



ai(/ci, cTi) + a2{k2, 0-2) -> asi^k^, 0-3) + a^i^k^, (T4) + Q{p, a) + Q{p, a) 



(3) 
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take the values cxj = ±1. The hehcity amphtudes for this process can be written as 



where the external quark or gluon currents are denoted by Jalai+2 ; ^^^^ the gluon propagator 
is given by 

DXS^d = ^ = 4 E ^M^^^ A,) , (5) 

for the conserved currents 

%',Jala^{ku h; (Xl, as) = q^'J^^aS^2, (^2, (T4) = 0. (6) 

We note in passing that the conditions in Eqn. [6] are satisfied not only for massless quark 
currents but also for the gluonic currents, in the processes (2b) and (2c) in the light-cone 
gauge [2] . The real gluon approximation to the amplitudes (Eqn. H]) is obtained by replacing 
the off-shell g*g* — )■ QQ amplitudes by the on-shell, and hence gauge-invariant gg — > QQ 
amplitudes 



e,.(gi,Ai)6,,(g2,A2)>l^;X^,„ (gi,g2,P,p;^,^) ^ (-^SJIfA, , (7) 



while keeping the four-momenta of the QQ system and the orientation of the colliding virtual 
gluon momenta in the QQ rest-frame. The full amplitudes for the process (3) are now 
approximated as 

A^o-CT ^ / 7-91 \Ai / jg2 \X2 ( kaQQ /o^ 

'o-io-3,(T2cr4 2„2 / ^ V'-^aia37o-i(T3V'-^a2a47o-20-4V-' '■gig2>>^iM ' \ ' 

where the incoming current amplitudes are given by 

(^a!W2)^:<..+2 = Ja.a,+2^^^, h+2; (^U ^i+2)e^(g^, A,)* (9) 

It has been shown in Ref. [2] that the helicity dependent phases of the above current am- 
plitudes reproduce azimuthal angle correlations among the jets while their magnitudes give 
the well-known DGLAP splitting functions when squared and summed over helicities. 
Following Ref. 0, we define the reduced current amplitudes j^^^-^^ as follows 

{'^alai+2)o'i<Ji+2 ^ ^3gl '^"^ Q iJ ulu ^+2 ' 

where Qi = \/—qf and the coupling factor ggl"'"'^^ is given by gsT^ for quark or anti-quark 
currents, and gsfabc for gluon currents. Here, gs = y/inag is the QCD coupling constant, 
while we suppress the colour indices in the helicity amplitudes of Eqn. H] for brevity. The 
reduced current amplitudes Jv/o-.+a have been computed and tabulated in Ref. [2] , to which we 
refer the reader for details. However, it might be worth recalling some of the salient features 
of these reduced current amplitudes. First of all, the reduced amplitudes for anti-quark 
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currents are same as those of quark currents. And since we neglect light quark masses, 
by chirality conservation, the helicity-fiip amplitudes will be zero for quark currents, i.e., 
J^Lfj = 0, while this is not the case for gluon current amplitudes. The existence of helicity- 
flip amplitudes in gluon currents give rise to certain special features in the amplitudes for 
QQ + 2 jets involving them. 

It is useful to look into the current amplitudes from the point of view of parton branch- 
ing [3], whereby the outgoing partons are emitted coUinearly. In the current amplitudes 
/ — )■ fg* or (yf — )■ gg*, if z is the energy fraction of the initial parton that is carried away by 
the off-shell gluon, then one can write the current amplitudes in terms of the energy fractions 
of the two t-channel gluons, zi, Z2, and the azimuthal angles of the outgoing partons (pi and 
02, in the frame in which the t-channel gluons are taken as propagating along the ±z-axis. 



2.1 Helicity amplitudes for gg — > QQ 

Let us first fix the reference frame in which we shall be working. As mentioned above, we 
take this frame to be the centre of mass frame of the two t-channel gluons producing the QQ 
pair, which is also the rest frame of the QQ pair. In this frame, the two gluons can be taken 
to be propagating in opposite directions along the z-ax.is with energy E each. Therefore, the 
scattering process, 

g{qi, Ai, ai) + g{q2, A2, 02) -> a, i) + Q{p, a, i) , (11) 

can be assumed to be taking place in the x — z plane. Here, ai, 02 denote the colour indices 
for the gluons and i, i denote the colour indices for the quark and the anti-quark. The four 
momenta of the gluons, Q and Q are given by 

gr = E(l, 0,0,1) 

g2^ = E(l, 0,0,-1) 

p^" = E{l,f3sme,0,Pcose) 

= E{l,-/3sme,0,-(3cose) (12) 



where /3 denotes the velocity of Q in the QQ rest frame and is given by /3 = ^ /l — , and 



QQ 

Mqq is the invariant mass of the QQ system. 

The helicity amplitudes for the process gg QQ can be expressed in the following form: 



( KA^Q \craii _ 2 



(13) 



One can further simplify Eqn. [I3]by noting the following relation 

MZ. = ^Zx,^P^ose . (14) 

Finally, the helicity amplitudes M'xlx^-, for different helicity combinations of the heavy quarks 
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and gluons, are given by 



MlZx = 2/3 sin 9{a\ + cos 9) /{I - (3^ cos^ 6) 

MZ = -2 v/r^(A + af3)/{l - (3^ cos^ 6) 
Mir = 2/3 a/1 -/3V sin^ e/{l - (3^ cos^ 6) 

= . (15) 

3 Azimuthal angle correlations: analytical results 

The squared matrix element (averaged over initial state and summed over final state colour 
(c) and spin (s)) can be obtained by using Eqns. |H1 [ISl ESI and the current amplitudes Jlla^^^ 
as tabulated in Ref. [2]. We find that for all of the different collision processes initiated by 
gg, qg and gg^ the squared amplitude can be cast in the following form: 

Vl/U- 16^gC (7/3 + 3/3^cos^g) ,2 • 2..g .2.2.. .4.^ 

+ 4/3^(1 - /32) sin^ ^{Fi cos(20i) + cos(202)} 

- 4(1 - /32)2F3 cos 2(01 - 02) - 4/3Sin^eF4 cos 2(01 + 02)] (16) 

where the colour factor C comes from the quark and gluon current amplitudes, averaged 
over initial colour, and is given by C = 1/36, 1/16 and 9/64 for gg, qg and gg sub-processes 
respectively. The co-efficients Fj's in the above expression are functions of the energy 
fractions of the initial partons carried away by the t-channel intermediate gluons. For nota- 
tional convenience, let us also define Zi = \ — zi. Also note that 0i and 02 are the azimuthal 
angles of the two outgoing partons (jets) in the frame defined in sub-section 12. 1[ 

The general expressions for Fq, which determines the total cross- sect ion, for the gg, qg 
and gg initiated sub-processes are as follows: 



l + zl^) (l + I^^) 

i+W^) (1+4+^^) 

^1 ^2^2 



j2) 

Z{Z\ Z2Z2 



The function Fi is found to be the same for the qg and gg initiated processes, and is given 
by, 

-l il + z^^) 

Hm] = ^ — -2 — 

^1 ^2 
^^(1 + 4 + ^) 



Zi Z2Z2 
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Similarly, the function F2 is the same for the qq and qg initial states, and is as follows: 



F r / 1 ^2 (1 + zi^] 



4 4 

Z2 {1 + Zf + Zf) 

= -2 ;2= ■ (19) 

Z2 Z-^Z\ 

The 2;j-dependent co-efficients of the terms proportional to cos 2(0i — ^2) and cos 2(0i + 02), 
namely F3 and F4, which will be seen to determine the azimuthal correlations in the threshold 
and the relativistic limits respectively, are found to be equal to each other (i.e., F3 = F4), 
and are also remarkably the same for all three sub-processes, 

F^Am/qg/gg] = • (20) 

We see that there is an enhancement in all the amplitudes in the zi ^ Q and ^2 limits, 
i.e., when the t-channel intermediate gluons become soft. Furthermore, in the amplitudes 
involving gluon currents (i.e., the qg and gg initiated processes), there is an additional 
singularity as — )■ 1 (i.e., — )■ 0). This happens because not only the intermediate gluons 
can become soft, but now an outgoing gluon jet can become soft as well. Since the form of 
F3 and F4 are identical for all of the gg, qg and gg collision processes, and these functions 
determine the angular distributions in the two kinematic limits to be considered next, the 
only difference that will arise is from Fq. We shall see that the ratio of F3 or F4 to Fq 
determines how strong the azimuthal correlations will be. 



3.1 Threshold and relativistic limits 



We now explore the two different kinematic limits in which one obtains characteristic strong 
azimuthal angle correlations that can be probed at the LHC The first one is when the QQ 
pair is produced near the threshold, i.e., the invariant mass Mq-q ^ 2mQ, and therefore, 
/3 — )■ 0. We can then express the squared matrix element in Eqn. [IB] in this simple form: 



1 2 /3^o, 



cicrs, (720-4 I 



Q1Q2 3 



cos 2(6^ — 



(21) 



From this equation, we can see that the differential distribution in (0i — ^2) will show a strong 
suppression at (0i — ^2) = 0, ivr, while it will show an enhancement at (0i — ^2) = ±vr/2 (we 
use the convention in which the azimuthal angles are in the interval [— 7r,7r]). Note that, in 
Eqn. |2T1 in the soft-limit for the t-channel gluons {zi — )■ 0), the combination AF^/Fq — )■ 1. In 
the next section we shall see that by studying the threshold production of top quark pair in 
ti +2 jets events at the LHC, we do observe these features in the corresponding differential 
distribution for (0i — 02)- We would like to point out that in the production of a CP-odd 
scalar boson in association with 2 jets one expects a similar distribution in (0i — 02) . 

The second kinematic region of interest is the relativistic limit in which the QQ pair 
invariant mass is very high, such that the masses of the heavy quarks can be ignored. In this 
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limit, as m <^ i?, /3 — 7- 1. Therefore, we can express the squared matrix element in Eqn. [TSl 
as follows: 



s.c 



\M 



(7(7 

(7ia3,o-20-4 I 



16^fC(7/3 + 3cos2e)(l + cos2 



QiQi 



sni' 



'e 



sin^ e 4F4 

+ C0S2 e Fq 



cos2( 



+ 



(22) 

From the above expression we see that in the relativistic limit, the differential distribution in 
the sum of azimuthal angles, (0i + (1)2), will show a strong suppression at (01 + 02) = 0, ±7r, 
while it will show an enhancement at (0i + ^2) = ±vr/2. In the next section, we shall study 
QQ + 2 jets {Q = b, t) processes at the LHC by demanding a high invariant mass for the QQ 
pair, and shall see that we are able to observe these features in the differential distribution 
for (01 + 02). 



4 Numerical results 

We now discuss the results for exact numerical evaluation of the differential distributions 
for (01 — 02) and (0i + 02) in the tt + 2 jets and bb + 2 jets processes. Here, in addition 
to the GF diagrams, all other diagrams at the tree level have also been taken into account. 
Using suitably devised cuts, we show that we are able to observe the angular correlations 
in different kinematic regions, as obtained using only the GF contribution in the previous 
section. We also compare the distributions and total cross-sections obtained by using the 
approximate matrix elements (presented in section |3l Eqns. [T6ll20|) with the exact ones. For 
exploring the threshold limit, we concentrate on the process tt-|-2 jets, with close to 2mt. 
For the relativistic limit, in which m E, we consider both the production of 66 + 2 jets 
and tt + 2 jets with a high invariant mass of the bottom or top pair. 

We present our results for LHC with a centre of mass energy of 8 TeV. The exact matrix 
elements were calculated with the help of MadGraph 5 ^4] at parton level. Subsequently, 
both the exact and approximate matrix elements were integrated over the relevant four- 
body phase-space using the BASES [5] Monte-Carlo integration package. We have used the 
CTEQ6L1 [B] parton distribution functions with LO running of a^, and as{Mz) = 0.13. For 
the factorization scale we choose the minimum value of transverse momenta demanded for 
the tagging jets {prmin ~ GeV). The strong coupling constant as is evaluated at the 
corresponding minimum transverse energy of the final state partons or heavy quarks, i.e., for 
the tt + 2 jets process, we take = aslnityasip'j-^^^y , while for bb + 2 jets it is taken to be 
«s = (^s{E^rnin)'^^s{pTmin)'^ ■ aualysis, -Ej^mm taken as 100 GeV (see Cut-3 below). 

The values of top and bottom quark masses used are 173 GeV and 4.7 GeV respectively [7]. 

To start with, we impose a set of selection cuts on the final-state partons (jets): 

Cut-1: 20 GeV < < 60 GeV, Ej > 250 GeV, \r]j\ < 5, ARj^j^ > 0.6 

Here p^, Ej and rjj denote the transverse momenta, energy and pseudorapidity of the jets 
respectively, while ARj^j^ {= ^ {rjji — rjj^y + (0ji — 0^2)^ j defines the separation between 
the jets in the pseudorapidity-azimuthal angle plane. We have subjected the tagging jets to 
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a slicing pj- cut, and have also demanded a high value for the minimum energy of the jets. 
Since we want to isolate the contribution of the GF diagrams to the full process, we also 
need to impose the so-called VBF cuts, which require the two tagging jets to reside in the 
opposite hemispheres of the detector and to be well separated in pseudorapidity. The VBF 
cuts imposed in our study are as follows: 

Cut-2: r]j^ > > r]j^, Ar]jj = r]j^ - r]j^ > 4. 

In addition to the above cuts on the tagging jets, we also impose a set of cuts on the t and b 
quarks. We demand a minimum transverse energy (E^ = ^Jm^ + for the h quarks, such 

that is always higher than the tagging jet Pt's. Also, the t and h quarks are required to 
be in the central region of the detector. These cuts are summarized below: 

Cut-3: E)^^ > 100 GeV, \v\b;b,t,t < 2-5 




Figure 2: Differential distribution of (0i — 02) (left panel) and (0i + 02) (right panel) for 
ti+2 jets in the three sub-processes initiated by qq, qg and gg and the sum of all sub-processes 
at the ^/s = 8 TeV LHC after Cuts 1, 2 and 3 as described in the text. The invariant mass 
of the ti pair has been demanded to be < 400 GeV. 



In Figure [21 we show the differential distributions obtained using the exact matrix el- 
ements for both the azimuthal angle difference of the two tagging jets, (0i — 02), and the 
sum of the azimuthal angles (0i -|- 02) in ti + 2 jets events, for the sub-processes initiated 
by qq, qg and gg. In addition, we show the sum of all sub-process contributions, too. In 
order to obtain the azimuthal angle distributions in the threshold limit, we need to ensure 
that the invariant mass of the tt system, Mu, should be very close to 2mt. Therefore, in 



8 



this figure, we have imposed an invariant mass cut of < 400 GeV. As explained in 
sub-section 3.1, we indeed see a distribution resembling the behaviour predicted by Eqn. [211 
In particular, we see the expected suppression at (0i — 02) = 0, ivr, and enhancement at 
(01 — 02) = ±vr/2. On the otherhand, the 0i + 02 distribution is observed to be rather flat, 
as expected in this region of Mu- This agreement between the analytic approximation and 
the exact numerical evaluation shows the usefulness of the VBF cuts in selecting the GF 
contribution to the total cross-section. In order to estimate the feasibility of studying the 
angular distributions in the tt + 2 jets process at the 8 TeV LHC, we show the cross-sections 
after cuts 1, 2 and 3, and various values of the cut in Table [T] As we can infer from 
this table, the total cross-section from all sub-processes taken together is considerable, such 
that with an expected integrated luminosity of ~ 20 fb^^, one can have sufficient number of 
events to study the azimuthal angle correlations at the 8 TeV run of LHC. 



Mtt cut 


aggtt (fb) CXggtt (fb) CTggtf (fb) ^44 + 2 j StS (fb) 


No Mtt cut 
Mtt < 400 GeV 
Mtt > 600 GeV 


655.80 1199.09 529.43 2384.32 
87.91 170.22 80.12 338.25 
195.62 334.58 137.42 667.62 



Table 1: Cross-sections of various sub-processes contributing to -|- 2 jets at 8 TeV LHC 
after Cuts 1, 2 and 3 and different values of the Mtt cut. We also show the total ti + 2 jets 
cross-section after these cuts. 




Figure 3: Normalized differential distribution of (0i — 02) (left panel) and (0i -|- 02) (right 
panel) for tt + 2 jets with different Mtt cuts, including all sub-processes, at the y/s = 8 TeV 
LHC, after Cuts 1, 2 and 3 as described in the text. The distributions obtained using the 
exact matrix elements are shown as histograms with solid lines, while the ones obtained 
using the approximate matrix elements are shown as curves with dotted lines. 
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In Figure [31 we show the normahzed differential distributions in (0i — and (01 + ^2) 
for the ti+2 jets process, where we present the sum of all sub-process contributions coming 
from qq, qg and gg. Here, we have shown the effect of variation in the Mu cut. While the 
01 — 02 distribution becomes more and more flat as we increase Mu, the 0i + 02 distribution 
is observed to resemble the prediction of Eqn. [22] for high Mjf. It is worth mentioning that 
in the distributions with no cuts on M^, one does retain some of the features of the two 
different kinematic limits. In this figure, we also compare the normalized distributions as 
obtained from the exact matrix elements and the approximate ones (the later being described 
by Eqns. [T6ll20|) . As can be observed, the distributions very closely agree with one another, 
thus confirming the validity of our approximation in correctly predicting the shapes of the 
azimuthal angle correlations. 



p-!p cut (GeV) 




'^Approx . / ^ Exact 






qqtt 


qgtt 


ggtt 


20 < < 60 


2.03 


2.02 


2.06 


10 < < 20 


1.08 


1.21 


1.39 



Table 2: The ratio of + 2 jets cross-section calculated using the approximate and exact 
matrix elements {<JApprox./o' Exact) for the various sub-processes, with different choices of the 
slicing cut for p^, at 8 TeV LHC All other cuts described in Cuts 1, 2 and 3 are kept fixed. 

We also show in Table [2]the ratio of the total cross-section obtained using the approximate 
and exact matrix elements. As we can see from this table, for the values of the kinematic 
cuts used so far in our analysis, although the on-shell gluon approximation predicts the 
normalized distributions rather accurately, it overestimates the total cross-section by about 
a factor of 2. Since our approximation approaches the exact matrix elements only in the limit 
when the t-channel gluons are collinear to the initial partons, we expect better agreement 
in the total cross-section when the p^'s are reduced further. This is what we observe in 
Table [21 where we see that by decreasing the p^s of the tagging jets, the total cross-sections 
also approach the exact ones. 

We now turn to the 66 + 2 jets process. Since the 6-quark mass is very small, we can 
easily obtain the characteristic distributions predicted for the relativistic limit by demanding 
a moderate cut on M^j. In Figure [H we show the differential distributions for both the 
azimuthal angle difference of the two tagging jets, (0i — 02), and the sum of the azimuthal 
angles (0i + 02) in 66 + 2 jets events, for the sub-processes initiated by qq, qg and gg. We also 
show the sum of all sub-process contributions separately. In this figure, we have imposed 
an invariant mass cut of Mj,^ > 400 GeV. In contrast to tt + 2 jets with < 400 GeV, 
the (01 — 02) distributions for 66 + 2 jets with the chosen value of M^j cut are found to be 
rather fiat, as expected. As far as the (0i + 02) distribution is concerned, as explained in 
sub-section 3.1, we indeed see a distribution as predicted by Eqn. [22] for the /3 — )• 1 limit. 
In particular, we see the expected suppression at (0i + 02) = 0, ±7r, and enhancement at 
(01 + 02) = ±7r/2. In order to estimate the feasibility of studying the angular distributions 
in 66 + 2 jets processes at the 8 TeV LHC, we present in Table [3] the cross-sections after 
Cuts 1, 2 and 3, and various values of the M^i cut. The total cross-sections for this process 
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Figure 4: Differential distribution of (0i — ^2) (left panel) and (0i + ^2) (right panel) for 
bb + 2 jets in the three sub-processes initiated by qq, qg and gg and the sum of all sub- 
processes at the \/i = 8 TeV LHC after Cuts 1, 2 and 3. The invariant mass of the bb pair 
has been demanded to be M^i; > 400 GeV. 



for different values of the invariant mass cut are seen to be considerably high, such that 
sufficient number of events for studying the angular correlations are expected at the 8 TeV 
LHC with ~ 20fb~^ luminosity, even after 6-tagging efficiencies are taken into account. 



Mft5 cut 


(^qqbb (fb) CTggbb (fb) aggbb ( ^3 ) cTfej+a jets (fb) 


No Mf,i cut 
Mf,i, > 400 GeV 
Mf,i > 600 GeV 


3748.27 7374.65 3593.10 14716.02 
1248.36 2333.83 1064.38 4646.57 
349.36 613.72 262.32 1225.40 



Table 3: Cross-sections of various sub-processes contributing to 66 -|- 2 jets at 8 TeV LHC 
after Cuts 1, 2 and 3 and different values of the M^i cut. We also show the total bb + 2 jets 
cross-section after these cuts separately. 

In Figure [5l we show the normalized differential distributions for the bb + 2 jets process, 
where we present the sum of all sub-process contributions coming from qq, qg and gg. Here, 
we have shown the effect of variation in the M^i cut. We also compare the normalized 
distributions obtained by using the approximate and the exact matrix elements, and they 
are found to agree to a good accuracy, as in the tt+2 jets case. We show in TablelHthe ratio of 
the total cross-section obtained by using the exact and approximate matrix elements for two 
different choices of the slicing pt cuts on the tagging jets. Here also, the approximate matrix 
element is found to overestimate the total cross-section, and the difference between the two 
is reduced once the tagging jet p^'s are demanded to be very small, thereby approaching the 
collinear limit for the t-channel gluons. 
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Figure 5: Normalized differential distribution of (0i — ^2) (left panel) and (0i + ^2) (right 
panel) for bb + 2 jets with different Mj,^ cuts, including all sub-processes, at the a/s = 8 TeV 
LHC, after Cuts 1, 2 and 3. The distributions obtained using the exact matrix elements are 
shown as histograms with solid lines, while the ones obtained using the approximate matrix 
elements are shown as curves with dotted lines. 



pip cut (GeV) 




Approx. / Exact 






qqbb 


qgbb 


ggbb 


20 < j4 < 60 


2.34 


2.31 


2.27 


10 < < 20 


2.07 


2.06 


2.07 



Table 4: The ratio of 66 + 2 jets cross-section calculated using the approximate and exact 
matrix elements {(TApprox./o-Exact) for the various sub-processes, with different choices of the 
slicing cut for p^, at 8 TeV LHC. All other cuts described in Cuts 1, 2 and 3 are kept fixed. 

5 Summary 

In this paper, we have studied the azimuthal angle correlations of jets produced in association 
with a top or bottom quark pair at the LHC. We have presented the helicity amplitudes for 
the process QQ+2 jets using only the gluon fusion diagrams and employing the on-shell gluon 
approximation to the matrix element. This gives us a general expression for the azimuthal 
angle distributions for processes initiated by qq, qg and gg initial states. Subsequently, we 
have explored this matrix element in two different kinematic limits. The first one is the 
threshold production of the QQ pair, in which the invariant mass Mqq is very near to 2mQ. 
In this limit we find a strong angular correlation in the distribution of the difference of 
azimuthal angles of the tagging jets, {(pi — (f)2)- The other kinematic limit considered is 
the relativistic limit, in which Mqq is very much higher than 2mQ. In this limit, while the 
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distribution in {(pi — is expected to be rather flat, there is a strong correlation in the 
sum of the azimuthal angles, {(pi + (p2)- We performed an exact numerical evaluation of 
these angular distributions in + 2 jets and 66 + 2 jets processes at the LHC, including all 
the diagrams at the tree level, and verified our analytical approximations using a suitably 
devised set of kinematic selection cuts. The VBF selection cuts are found to efficiently select 
the GF contribution to these processes. A comparison between the normalized differential 
distributions obtained by using the exact and the approximate matrix elements also shows 
a very good agreement between the two. The total cross-sections are overestimated by the 
approximate matrix elements by about a factor of 2, while they approach the exact cross- 
sections as the pr's of the tagging jets are demanded to be smaller, thereby achieving the 
coUinear limit for the t-channel intermediate gluons, in which our approximation approaches 
the exact matrix elements. 

We expect these azimuthal correlations calculated at the leading order to survive even 
after including higher order corrections, as has been previously demonstrated in the case of 
Higgs production with multiple hard jets [H [8]. One particular future direction that can 
be explored is to extend our parton level results by using a matrix elements matched with 
parton shower method. 

We would like to stress that the study of such angular correlations among jets produced 
with a top or bottom quark pair is not only interesting in itself, it has useful implications 
in the study of VBF processes for determining the spin and CP properties of new particles 
which can be discovered at the LHC. In particular, the experimental technique to measure 
such correlations between tagging jets can be established first by using these SM processes 
which have sufficiently high cross-sections. We expect that these studies can be performed 
with the already accumulated data at the 8 TeV run of the LHC. 
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